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INTRODUCCION 

Las ecuaciones diferenciales impulsivas se han comenzado a usar reciente- 
mente con exito para modelar el comportamiento de la biomasa de poblaciones 
individuales ante esquemas de captura pulsante masiva y peri'odos de veda. (ver 
P- ej [2]) 

En ciertos ambitos lacustres en que se permite de forma habitual la pesca 
deportiva tiene lugar un esquema distinto ya que la pesca se realiza de manera 
continua, y se incrementa en forma impulsiva la biomasa de una especie me- 
diante la siembra masiva de alevines de la misma. (Los Pescadores prefieren 
fuertemente a las truclias debido a la mayor resistencia que oponen) 

En los casos en que interactuan dos o mas especies competitivas bajo captura 
la dinamica se complica considerablemente. En rnuchos lagos patagonicos abier- 
tos a la pesca deportiva se produce dicha competencia entre especies autoctonas 
(existentes en la region incluso en el registro fosil), por un lado: las percas 
(Percichthys trucha -o trucha criolla, o perca de boca chica-, y Percichthys 
colhuapiensis -perca de boca grande-), y las especies exoticas: trucha arcofris 
(Oncorhynchus mykiss) y trucha rnarron (Salrno trutta), introducidas reciente- 
rnente por motivos comerciales. La practica de la siembra intensiva ocasional 
de alevines de truchas introduce variaciones pulsantes abruptas que facilmente 
pueden inducir una dinamica caotica(ver [3]). 

Como veremos, en determinadas condiciones las especies autoctonas llegan 
a veces a puntos cercanos a la extincion en determinado habitat. Si bien es 
posible reintroducirlas -aunque rara vez se hace-, no es lo ideal, ya que va en 
sentido contrario a la preservation de la diversidad genetica de dichas epecies. 

En otros casos el efecto de las siembras es menor al deseado, y la poblacion 
de truchas tiende a disminuir en promedio, a pesar de los impulsos. 

En este trabajo consideramos una estrategia para la decision de los mo- 
mentos de siembra teniendo en cuenta el valor del operador maximal de Hardy- 
Littlewood (a izquierda) para la rnasa de la poblacion a introducir; dicha estrate¬ 
gia puede presentarse como natural desde el punto de vista de quienes propulsan 
la siembra, y por otra parte podemos obtener en esta situation resultados que 
permiten controlar las magnitudes resultantes de las especies. 

Aqui consideramos ecosistemas relativamente cerrados (estanques, lagunas, 
lagos) donde puedan despreciarse las migraciones, en latitudes al sur del paralelo 
40 S -v.g., Lago Argentino, Laguna Verde, Lago Strobel, Lago Largo, Lago 
Futalaufquen, etc. En otras provincias argentinas: Cordoba, La Pampa, Buenos 
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Aires, suele realizarse tambien siembra de truchas pero allf el problema no se 
presenta debido a que en ellas hay peri'odos del ano en que las temperaturas 
relativamente altas, lo cual disminuye la cantidad de oxigeno disuelto en el 
agua. A partir de los 17/18° C las truchas tienen dificultades para digerir el 
alimento por carencia de oxigeno, y temperaturas superiores a 20° C resultan 
letales para las truchas. [4] 

Por simplicidad consideraremos el escenario donde interactuan 2 grupos de 
peces, ya sea porque solo hay una especie autoctona de interes y solo 1 exotica, 
o bien reuniendo por un lado las poblaciones autoctonas y por otro las foraneas: 

^, donde Ad es la biomasa del grupo originario (Pereas), y 

X 2 corresponde a los peces exoticos (Truchas). 

MODELOS IMPULSIVOS CON TIEMPOS DE IMPULSOS PREFIJADOS 

En las ecuaciones diferenciales impulsivas se consideran usualmente dos es- 
calas de tiempo, una continua donde la evolution es regida por una ecuacion 
diferencial, y otra discreta donde para una secuencia de instantes aislados la 
evolution es impulsiva: 

X'(t) = f(t,X(t)) si 
X(t + ) = g(k,X) si t = Tk 

Habitualmente los tiempos {• • • , r_ 2 , r_i, To, ti, T 2 , • • • } de impulso son pre- 
hjados (equiespaciados o no) y los impulsos X (r/ - ) dependen de k y solo del 
valor instantaneo de X (es decir X (r/ - )). 


X(t) = 


Xi ( t ) 
x 2 (t) 


MODELO IMPULSIVO CON TIEMPOS DE IMPULSOS 
DEPENDIENTES DE LOS VALORES HISTORICOS DE X 
Consideramos en cambio un modelo impulsivo en que los tiempos de impulso 
dependen de la imagen por X de los tiempos previos -asi corno tambien las 
magnitudes de los impulsos. 


/ X'(t) = f(t,X(t)) si t ^ T k 

\ X(t + ) = g(K,X ((— oo,Tfe])) si t = T k = infjt > r fc _ 1 : Q(X ((-00, t])) > X} 

MODELO CONCRETO CORRESPONDIENTE 

AL PROBLEMA DE PESCA.SIEMBRA 

En nuestro problema modelamos la situation en la escala continua por un 
modelo de competencia logi'stico en la biomasa total de ambas especies com- 
petidoras, donde K indica el valor de saturation, q un factor que indica la tasa 
de pesca, P un vector que indica el porcentaje no devuelto de las captura y 
C la matriz que indica las influencias rnutuas de las especies (la dependencia 
de X es solo a fin de la posible consideraciones de rnasas criticas para cada 
especie, debajo de las cuales cabe esperar que se extingan aunque sean no nulas 
-por ejemplo por escasas probabilidades de cruzamiento-). El segundo termino 
podrfa incorporarse al primero, modificando ( 7 , pero es preferible considerarlo 
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aparte para facilitar la interpretation. Los valores de C, q y P fueron tornados 
a partir de datos del comportamiento conocido de las especies, datos informa- 
dos de pesca y el relevamiento encuestado del aprovechamiento de los agentes 
respecto de las presas capturadas. 

X\t) = C (X(t)) *X(t)*(l- -q*(P,X (i)> 

Como las masas de las especies, y la contribution de las nrismas a la derivada 
son estimativas, y ademas pueden sufrir pequenas perturbaciones de diverso 
origen en realidad consideramos: 


X'(t) = c(x{t) s j *X(t) * 



- q*(p,X(t )) 


donde X = X + siendo 77 un vector aleatorio distribuido normalmente 
-cuya norma es relativamente pequena respecto de X. 

Luego, el modelo impulsivo clasico con impulsos equiespaciados sera: 


x'(t) = c(x(tfj *X(t) 


1 - 


X(t) 


K 


— q* ( P,X(t )) si t ± r k 


X(t + ) = g(k,X) si t = Tfc 


Como mencionamos, el exito de las politicas de siembra pulsante suele ser 
relativamente contingente, y frecuentemente ocurren situaciones de inestabilidad 
no deseadas, donde tiene lugar, a la larga, la prevalencia de una u otra especie, 
y esto para distintas tasas de siembra, y distintas frecuencias para la escala de 
tiernpo impulsiva. 

Esto se ha observado empiricamente en casos reales, y tambien puede notarse 
en modelos simulados, por ejemplo: 


MODELO CON TIEMPOS DE IMPULSOS DEPENDIENDO 
DE VALORES DE OPERADORES SOBRE X ((- 00 , t]) 

Teniendo en cuenta lo anterior hemos considerado una estrategia de siembra 
determinada por los valores de las biomasas. 


X'(t) = C * X{t) 


1 - 


X(t) 


K 


- q * ( P, A (t)) si t ± T k 


{ x ( t+ ) = o,r fc ])) si t = T k = inf{t > r fc _ 1 : Q(X ((- 00 , t})) > A} 

Es decir donde los impulsos (es decir las siembras) se realizan cuando una 
cantidad dependiendo de los valores previos de las masas de las especies supera 
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Figure 1: 


Figure 2: 


una determinada cota de tolerancia A. Y el valor del impulso depende tambien 
de valores actuales o previos, y posiblemente de la constante de saturation del 
habitat. 

Es conveniente considerar valores liistoricos, y no puntuales, ya que una 
medicion aislada no proporciona suficiente information para estimar la evolution 
del sistema. En la practica, dos sistemas similares que en un determinado 
momento registren iguales valores de X\ y X 2 en general no evolucionaran de 
la rnisina manera. 


Yendo al extremo, si en un instante dado se introducen en un estanque virgen 


arnbas especies en cierta proportion: X (to) 


X\ (t 0 ) 
X 2 (to) 


^, no se comportara 


igual que uno donde ambas especies conviven desde liace tiempo y en la que en 
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Figure 3: 


to la medicion tambien sea: X (to) = 


Entre otras causas esto se 


Xi (to) 

X 2 (to) 

debe a que la medicion releva involucra ejemplares adultos y se subestima la 
contribution de liuevas, alevines y ejemplares jovenes que mas adelante influiran 
en la biomasa. La medicion historica permite tener en cuenta dicha contribution 
ya que las cantidades de ejemplares inmaduros depende de ella. 

Un modelo propuesto, en linea con el interes de los lobbys de pesca, es 
considerar las diferencias entre la masa total de las percas y la de las truchas, 
y sembrar en el primer instante en que el mayor promedio de una ponderacion 
de las diferencias, en periodos previos supera A. 

Esto equivale a considerar para / (t) = a\Xi(t) — a 2 X 2 (t), el valor del oper- 
ador lateral (a izquierda de Hardy-Littlewood) M~ 


y el criterio para el impulso sera: 


t 

M~f (t) = sup J I f (s) ds > A 
h>0 n J 

t—h 

Es interesante el criterio, y tiene ventajas para el analisis debido a que es un 
operador de tipo debil (1,1). Por lo tanto si \f (t)\ tiende de rnanera razonable 
rapida a 0 (si las poblaciones tienden a equilibrar su diferencia en ciertos valores 
de coexistencia) esta acotada la rnedida del conjunto de tiempos en que M~ f ( t) 
supera A, en funcion de la norma 1 de /, es decir . 

En el caso del modelo de 2 especies competidoras -con siembra de alevines 
de una especie preferida (entre las 2 existentes)- nuestro criterio para decidir la 
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Figure 4: 

siembra depende de si el valor de un operador maximal (a izquierda, es decir en 
funcion de valores del pasado) supera o no un valor predeterminado A. 

Un criterio mas conservador apuntaria a realizar los impulsos solo cuando 
los promedios historicos desde el principio de las mediciones (y no el promedio 
mas favorable), a partir de un instante 0 arbitrario. esto equivale a considerar 
el operador de Hardy T : 

t 

Tf{i)=\ J f (s) ds 
0 

Tambien este es un operador acotado en L 1,0 °, y podrfamos considerar el 
criterio para el impulso: 

t 

Tf{t) = \ J f(s)ds> A 

o 

Esta claro, en cualquier caso, que M~ > T, aunque, corno dijimos, ambos 
operadores son acotados en L 1,00 , es decir que existe C > 0 : VA > 0 

\{t:M-\f\(t)>X}\<C\\f\\ 1 

Como M~ > T, la acotacion en L 1,oa de M~ implica la de T. Aqui 11 representa 
la rnedida de Lebesgue del conjunto. Pero tambien es posible probar la acotacion 
en un contexto muy general, con pesos, y con otras medidas. (vease por ejernlo 

[ 5 ]) 

Considerando / ( t ) = aiXi(t) — a 2 Af 2 (f), tal que f (t) —> 0 cuando t —» oo, 

es decir siendo — la relacion limite entre las dos especies competidoras si no 
a 2 
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Figure 5: 


hubiera perturbaciones, resulta que en este tipo de modelo (de Levins, ver [1]), 
resulta H/l^ < oo 

En la practica se pueden encuentrar objeciones a arnbos criterios: 
Consideramos, en arnbos casos los valores discretizados, con aj = f(j): 


M f (n) = sup (- a j ) 

n>k>l l n ~ k + 1 I 

1 n 

rm = -X> 

' j=1 

Vayamos ahora a la cuestion de que es lo que tiene que suceder para que 
el valor de uno de estos operadores supere un valor prefijado A en un instante 
n + 1, si liasta el instante n el valor del operador estaba por debajo de A: 

En el caso del operador maximal la cuestion depende fuertemente del nuevo 
valor a n +i- En efecto, supongamos 


A >M~f(n) 


1 

n — fco + 1 


n 

Y. a i 

j=k 0 


> 


1 

n — k +1 


n 

j=k 


(para cierto fco y para todo k). Ahora si M f (n+ 1) > A tiene cjue ser 

n +1 

M -/(n + 1 )= + i ^» j> A 

j = fc 

para cierto fc con n + 1 > fc > 1, luego 

n+1 

aj > {{n + 1) — fc + 1) A = A + (n — fc + 1)A 

j=k 
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Figure 6: 


, es decir 


, de donde: 


, por lo tanto: 


n 

«n+ 1 + E a,j > A + (n — k + 1)A 


n 

a n+ 1 — A > (n — k + 1)A — ^ Oj 

j—k 


cc n +1 — A ^ 

(n — k + 1) (n 


1 

ITT) 


n 

E a J 

j=k 


> A — M~ f (n) > 0 


, y entonces a n +i — A > 0, luego a n +i > A, que resulta una condicion necesaria 
para que M~f (n + 1) > A dado que A > M~ f (n). 


n +1 ' 

(ra+l)-fc+l E °-7 

, 3=k / 

que el caso k = n + 1 es menor que el caso general, es decir: 


Pero ademas, como M f (n + 1) = sup 

n+l>k>l 


resulta 


n+1 

\' 


\ v 7 v 7 -?=n+l / 


, de modo que si a„+i > A entonces A/ - / (n + 1) > a n +i > A. 

Luego M~ f (n + 1) > A. Es decir que si A > M~ f (n), entonces a n+ 1 > A 
es condicion necesaria y suficiente para que M~f(n+ 1) > A. 

Es decir, que en el caso discrete el criterio M~ f (n) > A, no difiere de 
observar la differentia ponderada puntual f (t) = ai-Xi(t) — CL 2 X 2 (t) e impulsar 



siempre que supere cierto valor A. La condition es entonces bastante laxa (no 
depende de n). 

Sin embargo uno quisiera tener en cuenta el comportamiento de la funcion en 
periodos relativamente extensos del pasado y no solo en el ultimo valor puntual. 
Consideremos por otra parte la situation para el operador de Hardy T. Si 

n 

\>Tf (n) = i a,j , para que 
j=i 


1 

n + 1 


n+1 

E< 

j=i 


= Tf (n+1) > A 


n+1 

necesitamos que aj > (n + 1) A, o sea 
i =i 


n 

Q>n+1 E aj > A + nX 

3=1 

n 

luego a n+ 1 — A > nA — a ? , entonces: 

j=i 


a„+i - A > n • (A 


1 

n 


E°j) 

j=i 


De aqui resulta que a ra _(_i — A > 0 aun es condition necesaria, pero en este 
caso ya no es suficiente. En terminos de la suficiencia debe ser 

a n +1 > A + n • (A - Tf(n )) 

que depende de n y de la diferencia entre Ay T/ (n). Es decir que el a„+i debe 
ser mayor que A en una cantidad igual a n veces la diferencia entre Ay T f (n). 
Por lo tanto es cada vez mas dificil al aumentar n, lo que a la larga representa 
una condition muy restrictiva. 


Decidimos entoncer usar como criterio que un promedio entre M~ y T supere 
A, es decir a ■ M~ f ( n ) + (1 — a) ■Tf(n ) > A, para balancear la puntualidad de 
M~ y la globalidad de T. 

Ademas, como 

{t : a ■ M~f (t) + (1 — a) ■ Tf (t) > A} C {t : a ■ M~ |/| (t) > A}u{f : (1 — a) ■ Tf(t) > A} 
se puede obtener la acotacion de a ■ M~ + {1 — a)-T a partir de las de M~ 

y t. 
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n+1 


Los valores de a n +\ y de cij claramente no son independientes, pero 

i — 1 

la covarianza tiende a 0 junto con n, de modo que si Tf (n — 1) es una fraction 

de A, es decir Tf (n — 1) = (1 — a) cA con 0 < c < para que el criterio se 

dispare es necesario, esencialenrente, que a„+i > acA. 

En las ensayos simulados, de hecho, con a = es decir con el criterio 

M~f(n)+Tf(n) . 

- 1 ----- > A, o sea, promediando sencillamente las condiciones ante- 


riormente consideradas obtuvimos un desempeno estable y bastante adecuado 
para distintos valores tasas de siembra. 

Cabe observar que, conro dijimos, cada uno de los operadores esta acotado 
en L 1,ao tambien el promedio lo esta, es decir existe C > 0 : VA > 0 


\{t : {M~ +T) I/I (t) > 2A}| < Cll/llj 

A continuation mostramos algunos ejenrplos de evoluciones, cuya estabili- 
dad era similar en varios centenares de ensayos (que incluian perturbaciones 
aleatorias normalmente distribuidas): 


Figure 7: 


CONCLUSION PROVISIONAL 

A modo de conclusion, los resultados de los ensayos sugieren que la es- 
trategia propuesta para la regulation de los poblaciones siguiendo el criterio 
de siembra determinado por los valores de M~ + T mantendria una evolution 
tendencialmente estable con distintas combinaciones de criterio de impulso y 
factor de siembra, y manteniendo acotada la magnitud total de las cantidades 
a introducir, permitiendo la pesca en los volumenes y frecuencias habituales y 
evitando las propensiones a la extincion. 

En principio, sobre los modelos numericos, los impulsos regidos por los op¬ 
eradores propuestos tenderlan en general a la estabilidad, ya que los impulsos 
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Figure 8: 


Figure 9: 

llevarian los estados de un ciclo en el espacio de fases a otro ciclo contenido 
en la capsula convexa del anterior. En la practica, hay que tener en cuenta 
que las especies pueden extinguirse en un determinado niclro ecologico en caso 
de que la cantidad de ejemplares caiga debajo de cierto umbra! critico y sera 
conveniente obtener datos de los resultados que se produzcan en el desempeno 
de estas estrategias en casos reales para sacar conclusiones definitivas. 
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Figure 10: 
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